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Abstract 

We construct a simple thermodynamic model to describe the melting of a supported metal 
nanoparticle with a spherically curved free surface both with and without surface melting. We 
use the model to investigate the results of recent molecular dynamics simulations, which suggest 
the melting temperature of a supported metal particle is the same as that of a free spherical 
particle with the same surface curvature. Our model shows that this is only the case when the 



q | contact angles of the supported solid and liquid particles are similar. This is also the case for the 



temperature at which surface melting begins. 
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INTRODUCTION 



Despite decades of study, the melting of nanoparticles continues to generate interest 
. In general, the melting temperature of spherical nanoparticles has been found 
to decrease in proportion to the surface area to volume ratio of the particle [5|, as the 
surface free energy of a molten droplet is less than that of the corresponding solid particle. 
Although free nanoparticle calorimetry has advanced considerably in recent years 

a a, 

most experimental determinations of the melti ng p oints of nanoparticles are conducted with 
supported particles (gold lead 9] and tin lQj . for example). The melting of supported 
nanoparticles is also important in carbon nanotube growth and other catalytic processes 
12], and for the stability of devices assembled from nanoparticles jisl. \h\. Thus it is of 
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interest to study the effect of the substrate on the melting point of a supported nanoparticle. 

Recent molecular dynamics simulations [ijj] of supported iron nanoparticles with a 
strongly interacting substrate found that the melting point of the particles was reduced 
in inverse proportion to the equilibrium surface curvature that results as they relax to wet 
the substrate. This statement also holds in the free particle limit since the curvature of a 
free spherical particle of radius a is proportional to its surface to volume ratio, 3/ a. Inter- 
estingly, the simulations in Ref. 15j found that the constant of proportionality between the 
decrease in melting point and the surface curvature did not depend on whether the particle 
was supported or free. In other words, the melting temperature of a supported particle that 
has a free surface with radius of curvature a, was found to be the same as that of a free 
spherical particle with the same surface curvature. The simulated nanoparticles in Ref. Ijj] 
also exhibited surface melting prior to complete melting. Surface melting is phenomena 
; nought to occur both on certain bulk metal surfaces 16j] and in certain metal nanoparticles 

13, U- 

In this paper we use a simple thermodynamic model to investigate the role of the sub- 



strate in both me 
the result in Ref. 



ting and surface melting of metal nanoparticles. Our model suggests that 



151 ]. that the relative decrease in melting point is proportional to the solid 



particle surface curvature, only holds when the contact angles of the supported solid and 
liquid particles with the substrate are similar. We also show that supported clusters will 
exhibit surface melting under certain circumstances, and that the surface melting tempera- 
ture in free and supported particles in clusters with same surface curvature is the same only 
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when the contact angles of the supported solid and liquid phases coincide. 



GEOMETRY OF SUPPORTED PARTICLES 

We start by considering a solid nanoparticle, initially spherical with radius a, that is 
placed on a flat substrate. We neglect the effect of faceting, curvature dependent surface 
energies and internal strains due to epitaxial mismatch with the substrate. Furthermore, we 
will assume that the particle has relaxed to its equilibrium geometry i.e. that the nanopar- 
ticle has relaxed to "wet" the substrate. Provided the nanoparticle is heated sufficiently 
slowly, the particle should relax to this geometry prior to melting. With these simplifying 
assumptions, the geometry of the relaxed particle can be approximated by a spherical cap, 
as shown in Fig. [Tj, with dimensions parameterised by either the cap height H, or radius of 
curvature of the free surface R, which minimizes the surface energy of the nanoparticle and 
substrate. 

The surface energy T of the system can be written as: 



where 7 S is the surface energy density of the free particle surface, 7& is the surface energy 
density of the substrate, 7 s b is the particle-substrate interfacial energy density and r& is the 
total energy of the bare substrate. We will assume that the density of the particle p s does 
not depend on the contact angle so that the volume of the supported particle remains the 
same as that of the free particle. 

Writing the volume of the particle as a function of H and R, it is then straightforward 
to show that T is minimized if H = — (A r y s b/ r y s ) R where Aj s b = 7b — 7 S — jsb- We note 
that Ajsb is often called the spreading parameter in the context of wetting phenomena: if 
A'jgb > then the particle will relax to fully wet the substrate. Here we are interested in the 
case where the particle does not fully wet the substrate (contact angles greater than zero) 
i.e. when Aj s b < and H/R = —A / ~f s b/ , ~f s > at equilibrium. In fact this minimum value 
of T can be written as 



where R* is the corresponding radius of curvature of the supported solid nanoparticle, given 



r = 2nRH ls + nH(2R - H) ( lsb - j b ) + T b 



(1) 




(2) 
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by 

(2 \ 2/3 ^ 
A^J (3 + A 7sb / 7s ) 1/3 ' (3) 
Thus T* and R* are the equilibrium surface energy and radius of curvature of the particle 

respectively. Note that the contact angle of the particle can range from to 180 degrees 
depending on the value of the spreading parameter A7 sfc . 

MELTING AND SURFACE MELTING 

In what follows we will assume that the density of the solid and liquid phases are identical 
i.e. p s = pi = p. We first consider the situation in which there is no surface melting. In this 
case, melting will occur at a temperature when the free energy of the solid particle wetting 
the substrate is equal to that of the corresponding liquid droplet wetting the substrate. 
If 7z is the surface energy density of the free liquid droplet and Rj is the corresponding 
equilibrium radius of curvature, then the difference in free energy between the solid and 
liquid is 

A M^H +3 (t-t)) (4) 

where f s (fi) is the bulk free energy density of the solid (liquid). Now, using fi — f s = 
pL (1 — T/T c ), where L is the latent heat of melting and T c is the bulk melting temperature, 
we find the melting temperature T m of the supported particle is given by: 

T m = T c (1 
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R* 



= ^W) + ^(l-|)T. (5) 

Thus if R* s = R] = R* then we recover the result of Ref. [id], namely that T rn = T^ ee (R*). 
In other words, if the contact angles of the solid and liquid droplets are equal, the melting 
temperature of the supported particle is the same as that of a free particle with an identical 
surface curvature, a = R*. However, if the curvature of the supported liquid particle is dif- 
ferent from that of the supported solid particle, it can be seen that the melting temperature 



will deviate from that found in Ref. 15] . 



Now we consider surface melting as illustrated in figure [T] which is thought to occur in 



many metals prior to melting 16]. We are interested in determining the onset of melting, 



when the solid particle is wet by a thin layer of melt (thickness 5) at the solid- vapor interface. 



We will assume that this melt forms a layer of uniform thickness with a geometry like that 
represented in figure [1] with 5 = R — r = H — h. The total free energy of the surface melted 
particle is then a function of 5: F(S) = V s (S)(f s — //) + V f\ + F(8) where V s ^ is the volume 
of the solid (liquid) and T is the thickness dependent surface energy. In particular 

T = 7T {2RH ll + r(2r - h) lsh + 5(2R - 5) llb + 2rh lsl {5)) 

where 7 s z(<5) = 7 S / + exp (—5/0 and £ is a correlation length that describes the thickness 



dependence of the interfacial energy in thin metallic liquid films 16| (in Pb, for example, £ 
has been measured to be ~ 0.6 nm [3]). As the surface melting proceeds, the curvature of 
the particle will relax to minimize the free energy i.e R* = R*(5) where R* minimizes the 
free energy F for a given S. 

In an isolated spherical nanoparticle of radius a, by minimizing the free energy F(S) with 
respect to 5 and setting 5 = 0, one can show that surface melting begins at a temperature, 
T s (a) given by 

T/ re » =T c (l-^ + 2 ^ s ~ t 1i A . (6) 

s K ' V p£ l P aL J 

provided A7 S ; > and a > £(7 S — jA/Ajsi (if a is less than this, full melting will precede 



surface melting i.e. T/ ree > T^ ree [17]] , and equation ([5]) will hold). 

For surface melting to occur in a supported solid nanoparticle with equilibrium curvature 
R s , a minimum in the free energy F(5) must appear at 5 = 0. It is straightforward to show 
that a minimum in F at 5 = occurs at the temperature T s : 



T s (R s ) = Ti ree (R 8 ) + 



P R S L I A 7 , 



fib 



T .free m , . H T c / COS fl B - COS flA 

- T ° {Rs) + -^L\ 1-cos^ J (?) 

where 9 S and 6 1 ; are the contact angles for solid particle and liquid particle respectively 
(defined via Young's relation 7^ cos^;) = % — Js(i)b)- Once again, if the contact angles of 
the solid and liquid droplets are equal, then the temperature at which surface melting occurs 
is identical to that of a free particle with the same surface curvature, R s i.e. T s = T/ ree (R s ). 
Further, if cos^ > cos 9i, so that the substrate favors contact with the solid over that with 
the liquid, the corresponding T s increases and vice versa. 

Complete melting will occur once the free energy of the surface melted particle, F(5), 
equals that of the corresponding liquid droplet, Fi i.e. at the temperature T m and liquid film 



thickness S m which satisfy F(S m ) = Fi. It is not possible to obtain an analytic expression 
for 5 m or T m , but numerical solutions to the resulting equations are shown in figure [2] as 
a function of R s for Pb particles. The figures clearly show the strong dependence of the 
melting temperature on the liquid droplet contact angle: a difference of ~ 10° in the molten 
particle contact angle can shift the melting point by ~ 100 K. Note that the melting point 
of a free particle with radius R s no longer coincides with that of a supported particle with 
radius of curvature R s when cos9 s = cos 9i, as in general the radius of curvature of the 
critical surface melted droplet will not be that of the solid particle (although the curves lie 
close to each other). 



CONCLUSION 



We conclude that the melting temperature (and surface melting temperature, if the par- 
ticle exhibits surface melting) of supported nanoparticles depends on the radius of curvature 
(or the contact angle) of both the supported solid and liquid particles. In general, we do 
not expect these curvatures to be the same: on a non-ideal solid substrate for example, 
epitaxial effects may favor one phase over the other. It is likely that the ideal substrate used 
in Ref 15J resulted in very similar solid and liquid particle contact angles. We have shown 
that is only in this "ideal" case that the melting temperature of free and supported particles 
with the same curvature is coincident, whether they exhibit surface melting or otherwise. 
Thus, results from free particle melting, where the curvature of the solid and liquid particles 
remain substantially the same, have only limited applicability to supported particle melting. 
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FIG. 1: The model for the geometry of a supported nanoparticle in equilibrium. We assume that 
the particle is a spherical cap of height H and radius of curvature R (left - the dashed lines simply 
illustrate the radius of curvature) . At the onset of surface melting, we assume that the geometry is 
close to that of the solid particle in its equilibrium geometry and that that the solid particle (radius 
of curvature r and height h) is initially wet by a molten layer of uniform thickness 5 = R—r = H — h 
(right). 




FIG. 2: The melting temperature T m and critical liquid film thickness 5 m for supported Pb clusters 
as a function of the radius of curvature R s of the relaxed solid particle for 7^ = 0.05, 0.10, 0.15 J 
m~ 2 and in the case where cos 9 S = cos Q\ (7^ ~ 0.13 J m~ 2 ). Also shown is the melting temperature 
of a free particle with radius R s . Other surface energies used are 7 S1 , = 0.61, ji v = 0.48, 7^/ = 0.05, 
7b = 0.25 and 7 s b = 0.1 J m~ 2 giving a contact angle of 75.8° for the solid supported cluster, and 
contact angles for the liquid droplets ranging from 78.0° to 65.4° respectively. Other parameters 
used were f = 0.63 nm, p = 10950 kg m~ 3 , L = 22930 J kg" 1 and T c = 600.65 K J 



